Scattering Theory of Nonlinear Thermoelectricity in Quantum Coherent Conductors 
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We construct a scattering theory of weakly nonlinear thermoelectric transport through sub-micron 
scale conductors. The theory incorporates the leading nonlinear contributions in temperature and 
voltage biases to the charge and heat currents. Because of the finite capacitances of sub-micron 
scale conducting circuits, fundamental conservation laws such as gauge invariance and current con- 
servation require special care to be preserved. We do this by extending the approach of Christen 
and Biittiker [Europhys. Lett. 35, 523 (1996)] to coupled charge and heat transport. In this way 
we write relations connecting nonlinear transport coefficients in a manner similar to Mott's relation 
between the linear thermopower and the linear conductance. We derive sum rules that nonlinear 
transport coefficients must satisfy to preserve gauge invariance and current conservation. We illus- 
trate our theory by calculating the efficiency of heat engines and the coefficient of performance of 
thermoelectric refrigerators based on quantum point contacts and resonant tunneling barriers. We 
identify in particular rectification effects that increase device performance. 

PACS numbers: 73.23.-b, 73.63.-b, 85.50.Fi 
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I. INTRODUCTION 

Today's energy challenges are calling for novel tech- 
nologies to produce, store and use energji These chal- 
lenges have among others renewed the interest in the 
old problem of thermoelectricity. Research in ther- 
moelectricity is predominantly directed towards finding 
new materials with better thermoelectric properties^—, 
however a recently proposed alternative approach has 
been to try and tailor low-dimensional micro- and nano- 
structured materials^, which may exhibit the sharp os- 
cillations in the electronic density of states necessary to 
trigger large thermoelectric effects^. With the ongoing 
trend towards circuit miniaturization in electronics, heat 
evacuation becomes a major problem and it is accord- 
ingly important to explore how heat currents can be con- 
trolled at sub-micron scales. Thus, there is ample mo- 
tivation for investigating thermoelectricity in micro- and 
nanodevices^. To name but few examples of recent in- 
vestigations, sub-micron electrically controlled heat rec- 
tifiers^, thermally driven electric rectifier o 10 ' 11 , thermo- 
electric refrigerators^— and heat converters^— have 
been proposed, and thermoelectric effects in molecular 
transport have been investigated^—. On a more funda- 
mental level, thermoelectric effects provide information 
on the nature of the underlying quantum fluid, for in- 
stance, breakdowns of the Wicdcmann-Franz law 2 ^ may 
reflect non Fermi liquid behaviors 2 ^—, quantum interfer- 
ence effect a 29 ' 30 or electric and heat transport mediated 
by both electron and hole quasiparticles in Andreev sys- 
tems2i~— . 

The field of thermoelectricity has until now over- 
looked the fact that thermoelectric devices often oper- 
ate in the nonlinear regime of transport (See for in- 
stance Refs. (Hull). 

Rectification effects occur in 
that regime^—, and their influence on thermoelectric- 
ity is not yet understood. There are only few works we 
know about which have discussed nonlinear thermoelec- 



tricity. Ref. [IjJ calculated the thermopower to quadratic 
order in temperature gradients and Refs. [42M44| focused 
on the nonlinear Peltier coefficient. These works con- 
sidered bulk systems with large capacitances, where ac- 
cordingly the build-up of local electrostatic potentials 
can be neglected. This neglect is no longer justified in 
small constrictions, where such potentials have to be self- 
consistently connected to the applied biases in order to 
preserve gauge invariance^. Refs. gH3 discussed non- 
linear thermoelectricity in such constrictions but did not 
introduce local potentials, thereby taking the risk of vio- 
lating gauge invariance. Refs. OUSil introduced sclf- 
consistcntly determined homogeneous Hartree potentials 
to preserve gauge invariance. 

In this manuscript we generalize the scattering theory 
of weakly nonlinear transport of Ref. [45[ to coupled elec- 
tric and heat transport. We go beyond Ref. [48[ which 
considered only electric currents in the presence of volt- 
age and temperature biases, by additionally treating heat 
currents. One of our motivations is to provide criteria 
for enhanced thermoelectric efficiency in the nonlinear 
regime. The standardly used figure of merit ZT is con- 
structed out of linear thermoelectric coefficients only 2 -, 
thus while any measure of thermodynamic efficiency in 
the linear regime can only be a function of ZT, it is 
expectable that ZT loses its predictive power in the non- 
linear regime. It is often the case, and we illustrate this 
below, that the maximal efficiency predicted by ZT is 
attained outside the linear regime. In this article we 
accordingly focus our attention on thermodynamic effi- 
ciencies and coefficients of performance as indicators of 
thermoelectric performance. In particular we investigate 
how they behave upon increasing sources of nonlinearities 
such as temperature and voltage biases and capacitive 
couplings. A key aspect of our theory is that in the non- 
linear regime of transport, the state of the system is self- 
consistcntly determined by the conditions under which 
transport proceeds. We will take this self-consistency 
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into account by introducing picccwisc constant local elec- 
trostatic potentials. We will go beyond Refs. [TilEIllIsj 
and consider more than one such potential. In the scat- 
tering theory, gauge invariance and current conservation 
are naturally preserved in the linear regime by the unitar- 
ity of the scattering matrix. In the nonlinear regime, they 
are reflected by sum rules connecting nonlinear tr ansp ort 
coefficients. Such sum rules were derived in Ref. [45j for 
purely electric transport and below we give a complete 
list of sum rules ensuring gauge invariance and current 
conservation for weakly nonlinear thermoelectric trans- 
port. 

The paper is organized as follows. In Section [Til we 
construct a scattering theory of weakly nonlinear ther- 
moelectric transport valid to quadratic order in temper- 
ature and voltage biases. We express linear and nonlin- 
ear transport coefficients in terms of the scattering ma- 
trix and the local potentials generated in the system by 
the biases. In Section IIII1 we discuss these local poten- 
tials and recall how they can be determined in practice. 
In Section IIVI we write down sum rules that have to be 
satisfied for the theory to preserve charge and current 
conservation as well as gauge invariance. In Section [V] 
we briefly discuss performance metrics for thermoelectric 
devices. In Section IVII we apply our theory to two par- 
ticular devices, a quantum point contact and a resonant 
tunneling double-barrier. We find in particular that the 
bias range of operation of these devices is such that our 
weakly nonlinear approximation gives a rather accurate 
description of the exact thermodynamic efficiency. Con- 
clusions are presented in Section IVIII 



II. SCATTERING APPROACH TO 
THERMOELECTRICITY 




FIG. 1: Sketch of a multi-terminal coherent conductor ca- 
pacitively coupled to an external gate and connected to three 
external Fermi liquid reservoirs. The electronic Fermi distri- 
butions of the reservoirs are depicted. In the nonlinear regime, 
an inhomogeneous electrostatic potential U(r) develops inside 
the conductor, as it is induced by the applied voltages and 
temperatures at the reservoirs. The arrows indicate our con- 
vention that electric and heat currents are positive when they 
flow into the conductor. 

In linear response, the scattering theory to thermoelec- 
tric transport was worked out by Butcher—. We follow 



his approach and consider heat and electric transport me- 
diated only by electron quasiparticles and extend it to 
the nonlinear regime. Considering a coherent conductor 
connected to i = 1, 2, ...M external reservoirs via reflec- 
tionless leads (see Fig. fTJ), the electric (Ii) and heat (Ji) 
currents through contact i are given by 

h = X E / d£ b ( £ ) A « t £ ' U W ' ( la ) 
J i = \Y,Jde{e-eV l ) f 3 (e) A l3 [e, U(f}] . (lb) 

We took the convention that both electric and heat cur- 
rents are positive when they enter the conductor from 
a reservoir. Thus, for example, the conductor is cooled 
if Ji < and heated otherwise. In Eqs. ([1]), e = 
E — Ep is the energy of an electron counted from the 
Fermi energy, fj(s) = {1 + exp [(e — eV^/ksTj]}^ 1 is 
the Fermi function in reservoir j, maintained at a voltage 
Vj and a temperature Tj, Aij [e, U(r)\ = Ni[e, U(r)]Sij — 
%j [e, U(r)] is given by the number iV, of transport chan- 
nel in lead i and the scattering matrix s via Tij = 
Tr[s|.Sij] where Sjj is the sub-block of s connecting ter- 
minals j to i and U(r) is the local potential in the con- 
ductor. In this manuscript, we take to be indepen- 
dent of energy and biases. The electrostatic potential in- 
duced within the conductor (see Fig.rjJ) is itself dependent 
on the voltages and temperatures in the reservoirs, i.e., 
U{r, {Vfc}, {T k }). Expanding the currents in Eqs. ([T]) to 
quadratic order in temperature differences, 9j = Tj — To, 
and voltage biases one obtains 

j 

+ i G i3kVjV k + B ljk 9 3 9 k + YijhVM , (2a) 

j-k 

Ji = ^ (TjjVj +5jj6j) 

3 

+ Y [TijkVjVk + Z, ,;.(),(), + * ijk YiO k ] . (2b) 

3,k 

The linear response coefficients were calculated in 
Ref. [5(| ■ They are given by 

Gn = \ I de(-dJ)A tJ (e), (3a) 

Bij = J^rJ dei-defleAijie), (3b) 

r« = f / de (-d e f) e Aijie) , (3c) 

Sij = Wol rf£ (-^/)£ 2 A v -( £ ), (3d) 

where f(e) = {1 + exp [e/fcgTo]} - . The temperature 
To 7^ is in principle arbitrary, however, better conver- 
gence of the expansion of Eqs. (HJ) is obtained when Tq 
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is taken as the temperature in one of the reservoirs or 
the average of the reservoir temperatures. In particular, 
the sum rules given below in Section [TV] guarantee that 
all currents Ii and Ji vanish at equilibrium, regardless 
of the choice of To. Eqs. (|3b[) and (|3c|) show that the 
reciprocity relation Tij = fly To is automatically satis- 
fied (see Ref. [5l| for a detailed discussion of Onsagcr 
relations in the linear regime). 

These linear response coefficients are determined by 
the equilibrium state of the system for which we set 
U^ eq \r) = 0. However, in the nonlinear regime, we must 
consider how the biases affect the local potential. We do 
that by introducing characteristic potential s] 45 ' 48 



U(r) = Y,lMr) Vk + z k {r) 8 k ] 



(4) 



which describe how U (r) depends on the biases to lead- 
ing order. Truncating the expansion of U (r) to linear 
order in the biases is sufficient to correctly capture cur- 
rents to quadratic order. Gauge invariance requires that 
Sfe u k(/) = l—i however there is no such relation involv- 
ing z k {r). The nonlinear response coefficients are then 

(5a) 

(5b) 

(5c) 



Gijk = J de (-d e f) a ijk (e) , 
2e 2 f 

Yijk = jTjrJ de (-d £ f) 7 4 jfc(e) , 
e 2 r 

Fijk = -r J de {-d B f) e a ljk {e) 



-G ik [ Si j - -5j,k 



(5d) 



Sijk = W 1 ] de { ~ dJ) £ Ajfc(£) + ^ Sj ' k ' (5e) 
2e f 

= jTjT J de {-def) e j ijk {s) 

-B tk {S itj - S jtk ) , (5f) 

where we defined 

6A 



aijk{e) = dr 



Pij k { e ) 



dr 



eSU{f) 
SAjjje) 
eSU{r) 



[2u k {r) - 5 jtk ] , (6a) 
[2eT z k {r) - sS j>k ] , (6b) 



li jk {e) = I dr I 6 -4^$ [eT Q z k {r) - eS j<k ] 



+ 



eSU{r) 
SA lk (e) 
e8U{r) 



euj{f) 



(6c) 



Several of the nonlinear response coefficients of Eqs. ([5]) 
differ primarily by a factor of e inside the energy inte- 
gral. Using a Sommerfeld expansion, we find that at low 
temperature the following relations between nonlinear re- 



sponse coefficients must hold, 

Fjjk + G ik (5i^ - $j,k/2) _ (7rfc B r ) 2 
Gij k 3e 
*ijk + B ik {5 id - S jtk ) _ (7rfc B To) 2 



E ijk - EijS jtk /{2T ) 



3e 



e ctln 



<9 £ ln G tjk , (7a) 
d £ \nY ijk , (7b) 

^ijSj, k 



2T n 



(7c) 



where all derivatives on the right-hand side are eval- 
uated at the Fermi energy, e = 0. These relations 
are the nonlinear counterpart to the Mott relation^ 5 - 
S = —{n 2 kgT /3e)d E InG between the thermopower S = 
— B12/G12 and the logarithmic derivative of the conduc- 
tance G = G12 in two-terminal conductors^. There is no 
counterpart to the Onsagcr reciprocity relations for non- 
linear coefficients, because, as can be seen from Eqs. ([5]), 
they all have different energy integrands. The absence 
of such relations may be attributed to the breaking of 
microscopic reversibility in the nonlinear regime, where 
local potentials are self-consistently related to transport, 
and where their dependence on temperature [the coeffi- 
cients z k in Eq. (QJ] is not necessarily related to their 
dependence on voltages [the coefficients u k in Eq. ((H)] . 



III. CHARACTERISTIC POTENTIALS 

Our main task is next to determine the characteristic 
potentials defined in Eq. (HJ. We briefly summarize how 
this can be don o 45 ' 48 . We consider first how additional 
charges are injected into the conductor as it is driven out 
of equilibrium. Biased reservoirs inject a net addition 
of charge, called the bare charge, Sq^. The presence 
of Sq^ induces an electrostatic potential shift within the 
conductor which, in its turn, generates a screening charge 
Sq^. The net injected charge is the sum of the two, 
5q{r) = 8q^ b \r) + Sq^{r). To leading order the bare 
charge can be expressed as a linear function of the biases, 



6q {b) {r) 



(8) 



where Dj{r) = e 2 / de(-dj) v?{r,e) and D^r) = 
e / de(—d E f) i/?(r, e) are charge injectivities (expressed 
in the units of capacitances as in Ref. [45|) determined 
by the energy dependent particle^ and entropion injec- 
tivities 



6s 



kj 



6s 



' fej eSU{r) eSU{r) 



s kj 



v e Ar,e) 



e 



i>j{r,i 



, (9a) 



(9b) 



where the trace runs over channel indices. Eq. ()9a|) gives 
the partial density of states at r corresponding to injec- 
tion from lead j, which is multiplied by a single-particle 
entropy factor e/T in Eq. (pb^ i 8 .. 
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Next, the screening charge 8q^ s \r) is related to U(r) 
via the Lindhard polarization function. The relation 
is in principle nonlocal, however the nonlocal correc- 
tions are quantum in essence and thus negligible in good 
metals with high enough electronic density. We thus 
make the quasiclassical approximation of a local Lind- 
hard function and write Sq^(r) = — H(r)U(r), with 
n(r) = Ej-D/Cf) = D (f)- From the definition of Dj 
we have D(r) = e 2 J de(~d e f) v{f, e) with the local den- 
sity of states v(r 1 e) = Ej i^(r, e). Note that at zero 
temperature, D(f) is the density of states (in units of ca- 
pacitance) at the Fermi energy. Putting all this together, 
the net total injected charge is, to leading order in the 
biases, 

Sq(r) = £ D)\r)V 3 + £ D](f)B d - D(r)U(r) . (10) 



The conductor may also be capacitively coupled to the 
reservoirs and external metallic gates. One way to take 
these couplings into account is to discretize the conduc- 
tor into different regions labelled r = 1, 2, each with a 
constant potential U r . The procedure brings about ad- 
ditional relations 



&Qr = E C rk(U r - Vfc) 



[11) 



with the sum running over all capacitively coupled com- 
ponents. This equation expresses changes in local charge 
distributions due to Coulomb interactions. The latter 
are represented here by geometric capacitances between 
different components of the system. 

Once the microstructure considered is consistently 
defined, in that all electric components (leads, gates, 
etc.) that are coupled to the conductor are included 
so that electric field lines do not leave the system, 
then the total charge is conserved by Gauss' theorem. 
Eqs. ([4]), pop and fTTj) then determine the character- 
istic potentials in Eqs. ([5]). The nonlinear transport 
coefficients are finally calculated by using the identity 
J dr 8u(jt\(...) = — e9s(—) to compute the functional 
derivatives in Eqs. ([6]). Two examples of thermoelec- 
tric microdevices will be discussed in Section EH which 
illustrate how the procedure is applied in practice by dis- 
cretizing the microstructures into regions of uniform local 
potential. 



IV. CONSERVATIONS AND SUM RULES 

Particle conservation is expressed mathematically by 
the unitarity of the scattering matrix, E, Aij = 
Y^j Aij = 0. This translates into the following sum rules 
for the linear transport coefficients, 

E G « = E B tj = r« = E H « = °- ( 12 ) 



In particular, two of these sum rules correspond to elec- 
tric current conservation, E, @ij = Ej Bij = while two 
others reflect gauge invariance, Ej Gij = Ej = 0. 
The unitarity of the scattering matrix further implies 

T,i a m = T,iPijk = EiTijfc = Ej (a»jfc + «ifcj) = 
y~]. jij). = 0, which results in the following sum rules 
between the nonlinear response coefficients, 



E G « fc 



= E B «* = E^* = ' 



(13a) 



E(g 



Gikj) = 



J 



3 



2_Yijh 



E*^ = ' 

j j 
E = ' 

i 

E r ^ fe ' 



G 



D 



{T ijk +T ikj ) = 0, (13b) 
(13c) 
(13d) 
(13e) 



Two of these sum rules were derived in Ref. [451 ] . 
Ei Gijk = 0, implied by current conservation and 
Ej (Gijk + Gikj) = 0, resulting from gauge invariance. 
All other sum rules are written here for the first time. 
Eqs. (|13a[) reflect current conservation while Eqs. (|13b[) 
and (|13c[) correspond to gauge invariance. The sum rules 

of Eq. (HiD imply Ei J% + I^Vi = 0, where /f n) is 
the linear response current through lead i. Thus in our 
approach energy is conserved to quadratic order in the 
biases. 

Finally we note the following, additional double-sum 



rules E 



jk a ijk 



j,k 



J2jk Ajfc = which imply 
Gijk = ^2 r «* = ' 



3,k 



Bijk — ^ijk 
j,k j,k 



(14a) 
(14b) 



While Eqs. l|14a|) are guaranteed to hold by Eqs. (I13bl) . 
Eqs. (|14b|) are new. Eqs. (JUJ), fI5c]) . (fT4aj) and (|14b|) 

ensure in particular that the electric and heat currents of 
Eqs. @ vanish at equilibrium. 



V. PERFORMANCE METRICS FOR 
THERMOELECTRICITY 

For two-terminal setups treated in linear response, 
Eqs. (P) reduce to 



G B 

r s 



V 

AT 



(15) 



where V and AT are the voltage and temperature dif- 
ferences between the reservoirs. Alternatively, this linear 
relation is often rewritten in terms of responses to an 
electric current and a temperature difference, 

j) = (tt p k) (at. 



(16) 
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Eqs. (|15j) and ([16)) define the linear response coefficients, 
the electric conductance, G, and resistance, R, the See- 
beck, S, and Peltier, ir p , coefficients as well as the ther- 
mal conductances S (at zero voltage, V = 0) and K (at 
zero current, 1 = 0). We can express the relationships 
between these linear transport coefficients as 



fl/G -B/G \ 

~ [t/g e-bt/g) 



(17) 



Dimensional analysis shows that a single dimensionlcss 
quantity can be constructed from the linear response co- 
efficients and the temperature. This is the figure of merit 
ZT = {GS 2 1 k)Tq. In the linear regime, thermodynamic 
efficiencies are then only functions of ZT and the tem- 
perature difference, for instance a heat engine working in 
the linear regime has maximal efficiency^ 



maxy 



IV 



AT\ ^JZT +1-1 



To J s/ZT + 1 + 1 



(18) 



It is important to keep in mind that this relation is valid 
only in the linear regime, and that nonlinear effects may 
invalidate it. Linear response may even predict a maxi- 
mal efficiency to occur outside its range of validity (at too 
large a bias voltage) . A description of thermodynamic ef- 
ficiency based only on ZT is therefore inappropriate in 
the nonlinear regime. However, the efficiency is generi- 
cally given by r\ = IV j Jhot , where Jhot is the heat current 
flowing from the hot reservoir. While Jhot = — Jcoid m 
linear response, Jhot 7^ — Jcoid in the nonlinear regime. 
Note that with this definition, it makes sense to substi- 
tute T) Jhot in Eq. (|18j) . in which case one obtains 
Carnot efficiency for ZT oo. 

The efficiency of thermoelectric refrigeration is mea- 
sured by the coefficient of performance (COP) which is 
the ratio of the heat extracted from the cold reservoir 
to the work done^, COP = J co \d/IV. This refrigeration 
"efficiency" is limited by COP < T co i d /AT. In the linear 
regime, its maximum is also given by ZT as 



COP m a X = max v 



J_ 

TV 



= (— 
\AT 



sJZT + 1-1 
y/ZT + 1 



1 

(19) 

where this time it makes sense to substitute To — ^ T co id. 
Note that the linear response scattering theory leaves the 
choice of To arbitrary, up to the fact that a well chosen 
To (close to both T co id and Tiot) minimizes nonlinear 
corrections. 

It is not at all guaranteed that ZT remains a valid 
figure of merit once nonlinear effects kick in. Below we 
discuss two examples of mesoscopic thermoelectric en- 
gines and compare their true efficiency to that predicted 
by linear response. 



VI. EXAMPLES 

A. Resonant Tunneling Barrier 

We illustrate our theory, first by calculating the ther- 
modynamic efficiency of thermoelectric devices based on 
a coherent single-level double-barrier system - a reso- 
nant tunneling barrier (RTB). Nonlinear electric trans- 
port through a RTB was investigated in Ref. [45j , where 
rectification effects were stressed in asymmetric systems 
with different tunneling probabilities through the left and 
right barriers. Electric transport generated by simultane- 
ously applied voltage and temperature biases was inves- 
tigated in Ref. [48j], focusing on the thermopower. This 
system is especially interesting in that the weakly nonlin- 
ear theory can be compared to fully nonlinear solutions^. 

r L r R 



ML ■ 



E r 



eU 



— Ef 

-—. - m 

Tr 



C 



FIG. 2: Sketch of a resonant tunneling barrier (RTB) formed 
between two Fermi liquid reservoirs by two tunnel barriers in 
series with transparencies Tl and Tr. They define a potential 
well with a single resonant level capacitively coupled to an 
external gate. 

The system we consider is sketched in Fig. [2] Two tun- 
nel barriers with tunnel rates Ti^r define a potential well 
with a single resonant level, at energy E r above the bot- 
tom of the well. Left and right, the well is coupled to two 
Fermi liquid reservoirs at temperatures Tl,r and electro- 
chemical potentials /iL,R- The well is capacitively coupled 
to an external gate with voltage V g which, together with 
the temperatures and electrochemical potentials in the 
reservoirs determine the potential U (assumed spatially 
homogeneous) in the well. This models a quantum dot 
in the Coulomb blockade regime as long as neither the 
temperatures nor the voltage biases exceed the single- 
particle level spacing in the dot. The equilibrium energy 
dependent particle injectivity at lead j =L, R is^ 



1 



7T (e + E F - E r 



r 2 



and the transmission probability is^ 

4r L r R 



T(e) = 



(e + Ep — E r ) 2 + r 2 



(20) 



(21) 
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where T = Tl + Tr is the total decay rate of the resonant 
level and we defined E F = + /Ur)/2. In the following 
we will set Ep = 0, as the zero of the energy scale and 
measure all other energies in units of T = const. 

Because we consider a single homogeneous potential in 
the well, U(r) = U, the injected charge loses its spatial 
dependence, Sq(r) —> 5q. Consequently, Eq. (fTU|) can be 
rewritten as 



(22) 



The well is capacitively coupled to an external gate, thus 
Eq. ([TTj) tells us that the additional charge in the poten- 
tial well is also given by 



Sq = C{U - V g ) . (23) 
and (1231) we obtain the characteristic 



Using Eqs. J22 
potentials as 



DY 



C 



C + D' 



C + D' 



C + D' 



(24) 



injectivitics 
i/^(e) and 



where the charge and entropic 
are given by Dj = e 2 J de{—d e f) 
D] = efds(-d e f) (e/T )^{e), with the particle 
injectivity i/| of Eq. (|20| . This fully determines the 
nonlinear coefficients in Eqs. (JS|) so that currents can be 
calculated using Eqs. (0). 

For the RTB the potential U can alternatively be de- 
termined nonperturbatively. Eq. (|22|) for the injected 
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FIG. 3: Electric (left) and heat (right) currents from a hot 
reservoir to a RTB with [E r — Ef)/Y — 1, /il — fJ.R = eV and 
Tr — Tl = 0.3 T . The base temperatures are k B T /T = 0.3 
(left) and ksTo/T — 1 (right). The different curves corre- 
spond to different tunnel barrier asymmetries (IT — Fr)/! 1 = 
0.95 (red), 0.5 (green), (blue), -0.5 (violet), -0.95 (cyan). 
Here, we switch off the capacitive coupling between the RTB 
and the external gate, C = 0. Fhe solid lines are the ap- 
proximate solutions in our weakly nonlinear approach and 
the dotted lines are the fully nonlinear solutions. The black 
dashed line is the heat current predicted from linear response. 
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FIG. 4: Left panel: Thermodynamic efficiency for a RTB- 
based heat engine with (E r — Ef)/T = 1, /il — /iR = eV, 
ksTo/r — 1, Tr — Tl = 0.8 To and different asymmetries 
(r L - r R )/r = 0.95 (red), 0.5 (green), (blue), -0.5 (vio- 
let), -0.95 (cyan). Right panel: COP for the same RTB as in 
the left panel with Tr — Tl = 0.1 To and working at larger 
bias, thus as a refrigerator. In both cases there is no capaci- 
tive coupling between the RTB and the external gate, C — 0. 
The solid lines are the weakly nonlinear approximate solu- 
tions and the dotted lines are the fully nonlinear solutions 
(see text). The dashed black line is the linear response solu- 
tion and the solid black line is the maximum efficiency/COP 
predicted from linear response. 



charge is replaced by the energy integral 

J2 u p j {e-eU)f j {e)-v{e)f{e) 



Sq = e / de 



(25) 



with v(e) = yT Vj(s). One determines U numerically 
from Eqs. (|2"3"|) and ([2~5"j) . The fully nonlinear currents 
are then obtained from Eqs. ([T]). The result is exact, 
up to the approximation of a single homogeneous local 
potential in the well. 

We first show in Fig. [3] the electric and heat currents 
for two RTB's. Different temperatures have been chosen 
to present / and J in order to visually optimize nonlin- 
ear effects. The curves are clearly nonlinear and display 
rectification in both the electric (left) and heat (right) 
currents. The rectification in the heat current is in gen- 
eral more pronounced, due to the additional voltage in 
the integrand of Eq. ()lb[) . Most remarkably, perhaps, 
we sec that the behavior obtained via the weakly non- 
linear approximation (solid lines) does not differ much 
from the fully nonlinear result (dotted lines), even for 
the strongly nonlinear heat current. The pronounced 
rectification effects will turn out to boost or hinder the 
efficiency of RTB-based devices depending on the asym- 
metry (r L - r R )/I\ 

We next investigate the efficiency of RTB-based ther- 
moelectric devices. Fig. 0] shows the thermodynamic ef- 
ficiency, 77 = IV/Jhot, °f a heat engine (left panel) and 
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the cocfficiccnt of performance, COP = J co \d/IV, of a 
refrigerator in the linear (dashed black line) and in the 
nonlinear approximation (colored solid and dotted lines) 
for different asymmetries (IY — r^J/P In linear re- 
sponse, I oc rLTpi/r and J oc TlTr, so that both r\ 
and COP are independent of (Pl — rpt)/r. For the cho- 
sen parameters, the maximal efficiency is at finite voltage 
bias, yet linear response gives a good approximation to 
r\ in the symmetric case, Tl = Tr. This is so, because 
in that case, Z? L ' e = and the screening potential 

given by Eq. (|2^f identically vanishes in the absence of 
gate coupling, C = 0. Introducing an asymmetry in the 
RTB changes the rules of the game, as rectification in 
the electric and heat currents kick in. As the asym- 
metry increases we can see a clear deviation from lin- 
ear response, which emphasizes the importance of treat- 
ing screening effects when characterizing thermoelectric 
performances of such microdevices. The efficiency 77 in- 
creases or decreases, depending on whether rectification 
play with or against the bias voltage. This depends on 
the sign of the bias voltage and the direction of rectifica- 
tion. From Fig. Owe see that at positive voltage, samples 
with large negative values for the asymmetry (IY — Tr)/P 
exhibit rectification effects where simultaneously the elec- 
tric current is enhanced while the heat current is reduced 
compared to their values in the linear regime. Both ef- 
fects conspire to increase the efficiency in that case (cyan 
curves in the left panel of Fig. The coefficient of 
performance, on the other hand, is always significantly 
smaller than predicted by linear response, and rectifica- 
tion plays in the opposite direction in the refrigerator 
case because COP = J co \d/IV while 7/ ~ IV/J\ wt . An- 
other remarkable feature is that both performance met- 
rics are well predicted by the weakly nonlinear theory 
(compare solid with dotted lines) with larger deviations, 
as expected, in the refrigerator case at larger biases. 

So far the capacitive coupling to the external gate has 
been set equal to zero, C = 0. Fig. [5] finally shows the 
dependence of 77 and COP on C. Efficiencies and max- 
imal efficiencies for a heat engine are shown in the left 
panels, while coefficients of performance and maximal co- 
efficients of performance for a refrigerator are on the right 
panels. For the chosen set of parameters, linear response 
overestimates 77 and COP when there is no coupling to 
the gate. In both panels c) and d), the capacitive cou- 
pling is seen to increase 77 and COP. Quite remarkably 
the efficiency ranges from well below the linear predic- 
tion for CT/e 2 = to well above it for CT/e 2 = 10" 5 , 
strikingly illustrating the dependence of 77 on nonlinear 
effects. The same calculation for RTB's with all other 
asymmetries reported in Fig. [4] show systematic increases 
in 77 above its value at CT/e 2 = as the capacitive cou- 
pling is turned on. The magnitude of the increase is in 
all cases similar to that in Fig. [SJ which we attribute to 
similar increases in U with C, regardless of the asymme- 
try. Panels a) and b) also show that better performances 
are obtained for a resonant level lying a distance of a 
few times T above or below Ep. We finally note that, as 
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FIG. 5: a) Maximal thermodynamic efficiency, 77 max = 
maxy?7, for a RTB heat engine and b) maximal coefficient 
of performance, COP ma x = maxvCOP, for a RTB refriger- 
ator as a function of the position of the resonant level, c) 
Efficiency and d) coefficient of performance vs bias voltage 
for RTB-based devices with (E r — Ef)/T = 1. In all cases, 
k B T /V = 1, (r L -T R )/r = 0.95 and the coupling to the gate 
voltage varies as CT/e 2 = (blue), 5 • 10 -6 (green), and 10 -5 
(red) with a gate voltage given by eV g /T = 3000. Panels a) 
and c) are for Tr — TL = 0.8 To and panels b) and d) are for 
Tr — Tl = 0.1 To. The solid colored lines are the approximate 
solutions in our weakly nonlinear approach and the dotted 
colored lines are the fully nonlinear numerical solutions. The 
dashed black line is the linear response solution and the solid 
black line is the maximum efficiency/COP predicted from lin- 
ear response. 



before and despite the rather large biases involved, the 
weakly nonlinear theory gives a rather faithful prediction 
of the true fully nonlinear behavior. 

From the data shown in this section we conclude that 
nonlinear effects can significantly increase the efficiency 
of RTB-based thermoelectric engines, first via rectifica- 
tion effects in both the electric and heat currents, and 
second via capacitive couplings which lift the resonant 
level. 



B. Quantum Point Contact 

Nonlinear electric transport through a quantum point 
contact (QPC) has been discussed theoretically in 



Ref. [45J and experimentally in Rcf. 53]. Linear ther- 
moelectric transport through a QPC has been discussed 
by Proetto^, assuming a saddle potential whose trans- 
mission spectrum is known exactly^—. Oscillations in 
the thermopower were in particular reported in Ref. [54| , 
which we reproduce in the inset to the right panel of 
Fig. [5] Curiously enough, the efficiency of QPC-based 
thermoelectric devices has not been calculated as of yet, 
even in the linear regime, and we briefly comment first 
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FIG. 6: Left panel: dimensionless ratio 3/GTLq for a QPC 
as defined in the text. This ratio is equal to one when the 
Wiedemann-Franz law is valid. Right panel: dimensionless 
figure of merit ZT. Left inset: conductance steps for the 
QPC. The arrows indicate the parameters corresponding to 
Figs. [8] and [9] Right inset: thermopower S of a QPC. Param- 
eters are CJ y /u x — 3 and T = 0.10 (blue), 0.25 (green), 0.50 
(red), 1.00 (violet). 



on the figure of merit for a QPC. We show in Fig. [5] 
data for the ratio S/GTLo with the Lorenz number 
L = 7r 2 (fc B /e) 2 /3 as well as ZT = {GS 2 /k)T . When 
the Wiedemann-Franz law holds, one has E/GTLo = 1—. 
This is the case when the QPC is well opened and trans- 
mits several conductance channels. The Wiedemann- 
Franz law breaks down when the QPC is closing, and 
the sharp increase displayed by S/GTLo m that regime 
suggests a reduction in ZT. The right panel of Fig. [S] 
however indicates an equally sharp increase of ZT reach- 
ing almost 10, which is understood once one realizes (see 
inset) that the increase in S/GTXo is accompanied by an 
increase of S^, with the net result that ZT increases as 
the QPC closes. This results in a large efficiency for an 
almost completely closed QPC operating as a heat en- 
gine. This however occurs when the power of the device 
is exponentially small. Some aspects of nonlinear ther- 
moelectric transport through a QPC have been addressed 
in Ref. [l4j]. We here extend these investigations. 

The geometry of the QPC is sketched in Fig. [7] Two 
lateral gates constrain the motion of electrons in a two- 
dimensional electron gas to a narrow channel. We take 
these two gates as symmetric and at the same potential 
and thus consider them together as a single gate. The 
width of the channel is tuned by the electric potential on 
the gates. We spatially split the QPC into a left and a 
right region, each with its own local potential. These two 
regions are capacitively coupled to one another (capaci- 
tance Co) and to the gate. For physically realistic param- 
eters, an asymmetry in the capacitive coupling of the two 
regions to the gate has only a small effect on transport, 
and we therefore consider a symmetric coupling with ca- 
pacitance C for both the left and right regions of the 




FIG. 7: Sketch of a quantum point contact (QPC) with the 
local potential in the conductor discretized into a left and a 
right region. The QPC is electrostatically defined by external 
gates. The capacitive couplings considered in our theoretical 
treatment are also indicated. 



QPC. We consider a geometrically symmetric QPC with 
only a few open transmission channels, which we model 
by a saddle potential^ 



eU(x,y) = eU 



1 



mujyy 2 



1 



2 2 



(26) 



with the electron mass m and the electric potential Uq 
at the center of the QPC. The transmission probability 
for each transmission channel is given 



T n (e) = [l + e- we » & 



with 

e«(e) 



E, 













-eU 







(27) 



ftu x , (28) 



for n = 0,1,2,... The density of states for each chan- 
nel, ^n(e), has been calculated in Ref. [58j], using a WKB 
approach. The result exhibits a nonphysical logarith- 
mic singularity close to conductance steps, but retains 
its validity away from them. We found that electric and 
heat currents lie on smooth curves except in the imme- 
diate vicinity of the steps. The results presented be- 
low accordingly use a numerical interpolation to inves- 
tigate the behavior of QPC-based devices close to con- 
ductance steps, whose exact form does not matter, given 
the very small range where the expression of Ref. [58| 
breaks down. Other than this interpolation, we use the 
density of states calculated in Ref. [58| . 

The QPC has been discretized into two separate re- 
gions, and we express the particle injectivity into region 
r = L, R from lead j =L, R as^ 



2[l-7;(e)]<y r 



(29) 



where the prefactor of 1/2 takes into account that only 
electrons moving toward the QPC undergo transmission 
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eV /hco eV / hco 

X X 

FIG. 8: Left panel: Efficiency of a QPC-based heat en- 
gine. Right panel: Coefficient of performance for a QPC- 
based refrigerator. In both instances the QPC is set at 
(Ef — eUo)/hw x = 1.5 corresponding to the red arrow in the 
inset of Fig. [6] Lu y /uj x = 3, eV g /hcj x — 3000, ksTo/huix = 0.3, 
Tr - T h = 0.1 T , eV = E F - fi R = ^ L - Co = and 
Chu) x /e 2 = (blue), 5 ■ 10" 5 (green), and 10" 4 (red). 




FIG. 9: Left panel: Efficiency of a QPC-based heat en- 
gine. Right panel: Coefficient of performance for a QPC- 
based refrigerator. In both instances the QPC is set at 
(EF — eUo)/fkdx = 0.5 corresponding to the black arrow in the 
inset of Fig. [6] ui y /cj x = 3, eV g /hoj x = 3000, ksTo/huJx = 0.3, 
Tr — Tl = 0.1 To, eV = Ef — MR — ML — /J-R, Co — and 
Chw x /e 2 = (blue), 2 • 10 -5 (green), and 4 • 10 -5 (red). 



or reflection. Eq. (fTU)) now becomes 

Sq r ~ J2 D rjVj + J2 D8 rA - DU r : 



(30) 



where r = L,R labels the left/right regions of the QPC. 
As in Eq. ([25)1 . the injected charge can alternatively be 
determined by 



6q r = C(U r - V g ) + C (U r - U P ), 



(31) 



where f = L(R) when r = R(L). Eqs. (|H]) and ^ 
determine the characteristic potentials as 



_ (C + Cp + D)D v r] + CpDfj C 



(32) 



z rg = , (33) 



(C + Cp + D)Drj + Co^ 
* rj *" (C + U)(C + 2C + U) ' * rg 



which, again, we use to calculate nonlinear transport co- 
efficients in Eqs. ([5]) and the currents in Eqs. ([2]). Note 
that the characteristic potentials now have two labels, 
one indicating the regions where they are measured, the 
other one the lead or gate from which they originate. 

We present our results for the efficiency of QPC-based 
devices. For all presented data, we set the chemical po- 
tential of the left terminal at Ef and bias the right reser- 
voir to lower electrochemical potential, = Bp — eV. 
This corresponds approximately to the situation consid- 
ered by Whitney^ (see his Fig. 2d), though he chose a 
Heaviside function for the QPC transmission. Another 
discrepancy between that work and ours is that we dis- 
cretize the geometry into two distinct spatial regions, 



which accordingly may carry two different local poten- 
tials. The data we present all have Co = 0, as we found 
that this capacitive coupling does not change the results 
significantly. Note that for the QPC, only the weakly 
nonlinear results are shown. 

In Fig. [5] we plot both the efficiency of a QPC-based 
heat engine (left panel) and the coefficient of performance 
of a QPC-based refrigerator (right panel) as a function 
of the applied bias voltage when the capacitive coupling 
C to the gate is varied. Both sets of data correspond 
to a first transmission channel close to half transmission, 
T — 1/2, corresponding to the red arrow in the inset 
of the left panel of Fig. |6l Fig. [8] shows that in the ab- 
sence of capacitive coupling, the maximal efficiency of 
the heat engine is slightly lower than one would predict 
from linear response, in agreement with Ref. [l4|. The 
trend is however reversed once the capacitive coupling C 
is turned on. Increasing C increases r\ and even extends 
the bias range of operation of the heat engine towards 
higher bias. Turning our attention next to the refrigera- 
tor case (right panel in Fig. [5J we note that the discrep- 
ancy between linear prediction and nonlinear calculation 
is stronger, essentially because the refrigerator works at 
higher bias voltage, therefore larger nonlinear effects may 
be expected. Increasing C increases the COP for refriger- 
ation just as it increased the power generation efficiency 
for the QPC when operating as a heat engine. In both 
instances we see that the nonlinear efficiency significantly 
exceeds the linear prediction at large enough capacitive 
coupling C. 

Fig. [S] shows a dramatic increase of ZT when the QPC 
is almost totally closed and its conductance is exponen- 
tially small. We investigate whether this region presents 
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opportunities for thermoelectricity, in particular whether 
the already large linear efficiency can be further increased 
by nonlinear effects. In Fig. |9] we plot data similar to 
those presented in Fig. [5] but obtained at smaller Ep, 
corresponding to the black arrow in the inset of the left 
panel of Fig. [6l We see that nonlinearities influence ther- 
modynamic performances only weakly in that regime, es- 
sentially because currents are exponentially small. 

VII. CONCLUSIONS 

We have presented a weakly nonlinear theory of cou- 
pled heat and electric transport in micro- and nanoscopic 
devices. On those scales it is of central importance that 
local electrostatic potentials are correctly taken care of, 
in order to preserve fundamental conservation laws such 
as gauge invariance and current conservation. We de- 
scribed how this can be done and listed sum rules en- 
suring that these conservation laws are preserved. Most 
works on thermoelectricity use the linear figure of merit 
ZT as a metric for thermodynamic performance - the 
larger ZT, the better the efficiency of a heat engine 
and the coefficient of performance of a refrigerator. We 



pointed out that the predictive power of ZT is no longer 
guaranteed when nonlinear effects kick in, and directly 
investigated the thermodynamic efficiency and coefficient 
of performance, including nonlinear contributions to the 
heat and electric currents. We investigated devices based 
on two microconstrictions, a quantum point contact and 
a resonant tunneling barrier. In both instances we found 
situations where the true efficiency significantly differs 
from that predicted by linear response, as well as cases 
where thermodynamic efficiencies strongly depend on the 
strength of these nonlinearities - not only on the magni- 
tude of the applied temperature and voltage biases, but 
also on the capacitive coupling to external electrostati- 
cally charge regions. Most importantly, we found that 
the weakly nonlinear theory faithfully captures the true 
nonlinear behavior of micro- and nanoscopic devices. We 
believe this opens the way to more general investigations 
of thermoelectricity in quantum microdevices which until 
now have often been limited to numerical investigations. 
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